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Abstract. 

Terrestrial planetary systems may exist around nearby stars as the Earth- 
sized counterparts to the many giant planets already discovered within the 
solar neighborhood. In this chapter we first discuss the numerous techniques 
which have been suggested to search for extrasolar terrestrial planets. We 
then focus on the expected results from that technique in which an orbiting 
telescope or interferometer is used to obtain a visible or infrared spectrum 
of a planet, without contamination from the parent star. We show examples 
of such spectra for selected cases: the present Earth, the Neoproterozoic 
(snowball) Earth, a methane-rich Earth, and the present Mars and Venus. 
We conclude by discussing the implications of such spectra for the detection 
of life on an extrasolar terrestrial planet. 

To appear in "Atmospheres in the Solar System: Comparative Aeronomy", 
edited by M. Mendilo, A. Nagy, and H.J. Waite, AGU Geophysical Monograph 
130, 369-380, 2002. 



1. KNOWN EXTRASOLAR 
SYSTEMS 

The first planet orbiting a solar-type star beyond 
the solar system was announced in 1995. Since then, 
as of October 2001, 66 planets have been found and 
confirmed, orbiting 58 stars with a median distance 
from the sun of about 28 pc. The search database 
contains roughly 1200 stars. 

Current estimates of the frequency of massive plan- 
ets range from about 3-5% to 6-7% [J. Schneider, and 
G. Marcy, resp., personal communication]. The num- 
ber of detected planets is growing monthly, as ob- 
serving techniques are refined, and as the time base 
of the record increases, allowing longer-period planet 
signatures to be extracted from the radial velocity se- 
quences. 

The definition of a planet can be a controversial is- 
sue, as is evidenced by the recent debate on whether 
Pluto should be classified as a planet, or a trans- 
Neptunian or Kuiper-belt object. However in the case 
of extrasolar planets the debate centers not on the 
question of the low-mass end of the scale, as for Pluto, 
but on the high-mass end, which is the realm of the 
planets discovered to date. The question is, should 



these objects be classified as planets, brown dwarfs, 
or small stars? 

Here we adopt the recommendation of Oppen- 
heimer et al. [2000] for objects of solar metallicity: 
the minimum mass for a brown dwarf is about 13 Mj 
(where Mj is the mass of Jupiter) or 0.013 M su „, 
sufficient to allow deuterium burning; the minimum 
mass for a main sequence star is about 78 Mj or 
0.075 M su „, sufficient to allow hydrogen burning. 

Pulsar planets, known since 1992 from variations 
in pulse arrival times from pulsars, are kept in a sepa- 
rate category from exoplanets around main-sequence 
stars. These bodies may have formed during the ex- 
plosion that created the neutron-star pulsar. 

At least two websites collect current information 
on exoplanets. One site is "The Search for Extraso- 
lar Planets" [Marcy, 2001] originating at UC Berke- 
ley. The other, larger, site is "The Extrasolar Plan- 
ets Encyclopaedia" [Schneider, 2001] originating at 
Paris Observatory, Meudon. Both sites feature news 
items, discussions, tutorials, papers, bibliographies, 
and comprehensive lists of exoplanets, their proper- 
ties, and their parent stars. Both sites are expertly 
edited and authoritative. 

The present paper is oriented toward extrasolar 
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terrestrial planets, i.e., those in the mass and tem- 
perature range of Venus, Earth, and Mars. However, 
since all known exoplanets are in the gas giant mass 
range, it is appropriate to ask about the prospects for 
finding terrestrial-mass exoplanets. 

Astrometry (cf. later section) tells us the semi- 
major axis, eccentricity, and mass of a planet. Each 
of these is a clue to a part of the history of the planet, 
e.g., where and how it was formed, whether it has 
migrated since then, and whether its orbit has been 
perturbed by other bodies. However to learn from 
this data if terrestrial-sized planets might be present 
we must fall back on inference based on the observed 
frequency distribution of mass, as follows. 

If we examine a plot of exoplanet mass versus semi- 
major axis, we see that the discoveries tend to pop- 
ulate the entire region between the extremes of: (a) 
maximum mass set by definition at about 13 Mj; (b) 
minimum mass set by the radial velocity detection 
method at about M/Mj = 0.035-ui? 1/2 where v (m/s) 
is the minimum detectable orbital velocity amplitude 
and R (AU) is the orbital radius; (c) maximum ra- 
dius set by a total observing time of about 5 years, to 
get a full orbit; and (d) minimum radius set by the 
smallest observed orbit at about 0.01 AU. If we as- 
sume that this entire region, on a (log i?, log M) plot, 
might be uniformly populated, then we can calculate 
the expected number of planets as a function of mass. 

We show in Figure |l| a histogram of the detected 
masses, and two "expected" distributions correspond- 
ing to velocity amplitudes of 2 and 20 m/s, the former 
being a nominal goal of the present radial velocity 
searches, and the latter being an estimate of the cur- 
rent level of confident detection. Although Figure [l]is 
undoubtedly overly simple, we may nevertheless draw 
two conclusions. (1) At the high-mass end, where the 
present searches are certainly relatively complete and 
unbiased, the observed distribution drops off rapidly 
from about 1 Mj to 13 Mj, suggesting that the exo- 
planet population really does exist as a separate entity 
from any brown dwarf population (not shown) which 
might fill in the 13-75 Mj range. (2) At the low- 
mass end, the observed population seems to be ap- 
proximately fit by the 20 m/s curve, suggesting that 
the observations are observationally limited, and that 
there is no evidence for a fall-off at low masses. Thus 
at present we may well be sampling only the very top 
end of the exoplanet mass distribution curve, and it 
may well be that the distribution function continues 
all the way down to the M-Earth range. On this ba- 
sis we optimistically look forward to someday finding 
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Figure 1. Histogram of discovered exoplanets, rang- 
ing from about 0.16 to 13 Mj. Theoretical curves are 
shown for the cases where exoplanets have a uniform 
density distribution in the (logM, logi?) plane, in the 
range logAf = —3 to +1.11, and logi? = —2 to +1, 
and limited by the sensitivity of radial velocity mea- 
surements to orbital velocities of 2 and 20 m/s, as 
indicated. The observed distribution appears to be 
bounded by measurement accuracy on the low-mass 
side, and by a lack of exoplanets on the high-mass 
side. 

terrestrial exoplanets. 

2. DETECTION METHODS: 
GRAVITY AND AERONOMY 

There is a surprising number of proposed tech- 
niques to detect or characterize exoplanets. About 
one-half of the techniques may be classified as being 
gravitational in nature, and the other half as being 
aeronomic (or photonic). The gravitational methods 
are (1) radial velocity, (2) astrometry, (3) transits, 
(4) pulsar timing, (5) gravitational lensing, and (6) 
disk shaping. The aeronomic methods are (7) visible- 
infrared shift, (8) reflected light, (9) transmitted light, 
(10) auroral emission, (11) radio emission, (12) an- 
thropogenic transmission, (13) coronagraphic imag- 
ing, and (14) interferometric imaging. 
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2.1. Radial velocity 

As the planet and star orbit their common center 
of mass, the velocity vector of the star projected along 
the observer's line of sight is proportional to Msinz 
where M is the planet mass and i is the inclination 
of the orbit plane to that of the sky. Doppler shifts 
in stellar spectra have been measured to an accuracy 
of about 3 m/s; 10 m/s is common, and 1 m/s may 
be the ultimate limit of this technique [Marcy and 
Butler, 2000; Santos et at, 2000]. By comparison, 
the solar velocity due to Jupiter is about 3 m/s, and 
that due to Earth is about 0.01 m/s. 

2.2. Astro metry 

The projection of stellar orbital motion onto the 
plane of the sky produces an astrometric shift which 
is measured with respect to a grid of nearby refer- 
ence stars. At 10 pc, Jupiter would move the Sun 
by about 100 /xas (1 ^as = 10~ 6 arc-sec), and Earth 
would move the Sun by about 0.3 /xas. The Hip- 
parcos satellite had an accuracy of about 500 /ms, 
so could not quite detect Jupiter-sized planets, but 
planned missions such as FAME (50 /uas) [Horner et 
al, 2001] (2004 launch planned), GAIA (2-10 ^as) 
[Perryman et at, 2001], (2012 launch planned), and 
SIM (1 /ias) [Danner and Unwin, 1999] (2009 launch 
planned), have a good chance of detecting masses in 
the sub- Jupiter range, and almost down to the terres- 
trial limit. 

2.3. Transits 

If the planet's orbital plane is seen nearly edge-on, 
a partial eclipse of the star by the planet may oc- 
cur. For the one example known to date, HD209458, 
[Charbonneau et al., 2000], precise photometry has 
allowed us to infer the stellar limb darkening, the 
planet radius, the orbital inclination, and therefore 
the planet's mass. Further observation may lead to 
a transmitted light measurement, as discussed below. 
Dedicated, staring telescopic searches [e.g., Borucki 
et al., 2001] may detect more examples of this rare 
transit phenomenon. 

2.4. Pulsar Timing 

The clock-like constancy of pulsar spin rates means 
that the time delay produced by line of sight dis- 
placement due to an orbital companion can be in- 
terpreted in terms of the mass and orbit of the com- 
panion [e.g., Konacki et at, 2000]. Two such pul- 
sars are known, PSR 1257+12 (3 planets, possibly 4), 



and PSR B1620-26 (1 planet). Interestingly the tech- 
nique is sufficiently sensitive that we know that some 
of these pulsar planets are in the few MEarth range, 
and one may be in the Mpi uto range. 

2.5. Gravitational Lensing 

A star in our Galaxy can cause a distant light 
source (a background galaxy) to apparently brighten 
for several days by gravitationally deflecting and lens- 
ing the distant galaxy's light as the star happens to 
pass in front of the galaxy. If the star hosts a plan- 
etary companion, then the planet can cause a sec- 
ondary brightening [Gaudi and Gould, 1997]. The ef- 
fect would be transient, but statistics of many events 
would give information on the Galactic incidence of 
planets. 

2.6. Disk Shaping 

For a star with a debris disk, the presence of a 
planet with non-zero eccentricity or relative orbit in- 
clination will cause the disk to become eccentric or 
warped, respectively; a planet inside the dust or de- 
bris disk can also generate resonant trapping (clumps) 
and clearing of central holes [Wyatt et al, 1999; Liou 
and Zook, 1999]. 

2.7. Visible-Infrared Shift 

Precise astrometry of an unresolved star-planet 
system in the visible and infrared wavelength regions 
simultaneously would show a spatial shift because the 
star light dominates the planet's light less in the in- 
frared. This is a detection mode planned for the Keck 
Interferometer and Very Large Telescope Interferom- 
eter [Akeson et al, 2000; Lopez et al., 2000]. 

2.8. Reflected Light 

For a large, close-in planet, sufficient star light may 
be reflected from the planet that it could be directly 
detected as an additional component of star light, 
even if the system is not spatially resolved [Charbon- 
neau et al, 1999]. 

2.9. Transmitted Light 

If a planet eclipses the parent star, then the 
planet's transparent upper atmosphere will trans- 
mit a portion of the starlight, with a superposed 
planetary absorption spectrum which could produce 
measurable features such as He 1.0830 /xm or Na 
0.589 /im [Brown et al, 2001]. 
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2.10. Auroral Emission 

Auroral activity such as seen on Earth or Jupiter 
(e.g., O 0.5577 /im, Ha) is a non-thermal, potentially 
useful indicator of a planet [e.g., Waite et ai, 2001], 
but the flux rate may be so low as to not be com- 
petitive with other techniques. Detection of oxygen 
emission is discussed more fully in Chapter VI. 2. 

2.11. Radio Emission 

Decametric radio wavelength radiation from elec- 
trons in the magnetic field around Jupiter and Io sug- 
gests that planets might be identified by this non- 
thermal radio signature [Bastian et ai, 2000], but low 
flux might limit this method. 

2.12. Anthropogenic Transmissions 

(a) It has been pointed out that a directed visible 
laser beam from Earth could be made to outshine the 
sun for the duration of a pulse, and would therefore 
be easily visible with a modest telescope at interstel- 
lar distances [Howard and Horowitz, 2001]. In the 
hope of detecting such pulses trial experiments have 
been started, (b) Pulsed radio transmissions contain- 
ing coded messages are likewise detectable at great 
distance, and have been the basis of several searches 
[e.g., Horowitz and Sagan, 1993; Leigh and Horowitz, 
1999]. 

2.13. Coronagraphic Imaging 

Direct visible-wavelength detection of the analog 
of the solar system's gas giant outer planets could 
be achieved using existing general-purpose telescopes, 
such as the Hubble Space Telescope or the planned 
Next Generation Space Telescope, if the residual op- 
tical imperfections in either telescope were to be 
corrected by adaptive optics. Dedicated corona- 
graphic telescopes with shaped or shaded pupils have 
also been proposed [Nisenson and Papaliolios, 2001; 
Spergel, 2001]. To detect an Earth at 10 pc and 0.5 
/jm wavelength requires 10~ 10 starlight rejection at 
0.1 arcsec separation. 

2.14. Interferometric Imaging 

An infrared- wavelength imaging interferometer has 
been proposed as a NASA mission [Beichman et at, 
1999], or European Darwin mission [Fridlund, 2000]. 
To detect an Earth at 10 pc and 10 fim wavelength 
requires 10 -7 starlight rejection at 0.1 arcsec separa- 
tion. 



3. WHY SPECTRA OF 
TERRESTRIAL PLANETS? 

Of all the techniques mentioned to detect or char- 
acterize an extrasolar terrestrial planet, we have cho- 
sen to focus on just two: coronagraphic imaging and 
interferometric imaging. The reason for this is that 
these seem to give us the best chance to determine 
the atmospheric constituents of the planet, by direct 
observation of the reflected or emitted light from the 
atmosphere and surface. 

The current thumbnail picture of planetary forma- 
tion is as follows. A massive molecular cloud is some- 
how triggered to collapse; a star is formed; the star is 
surrounded by a remnant gas and dust cloud; rocky, 
metallic, and, in the outer, colder parts of the cloud, 
icy grains, condense and agglomerate in the surround- 
ing cloud; and the agglomerations cascade to larger 
sizes to form planetary cores. Then in the outer part 
of the cloud, where plenty of gas is available, the gas 
continues to collapse around the cores, and Jupiter- 
like (gas giant) planets form. In the inner part of 
the cloud where it is hotter and less gas is available, 
only rocky planets form. These rocky, terrestrial-type 
planets can have abundant liquid water or ice, and 
relatively thin atmospheric envelopes, both generated 
by outgassing of the rocky material and possibly also 
by infalling comets. It is on these planets, with their 
solid surfaces where water can accumulate and chemi- 
cal reactions occur, and their thin atmospheres where 
sunlight can penetrate and be used for driving chemi- 
cal reactions, that we speculate that life probably has 
its best chance to develop. (Excellent references can 
be found in Chyba et al. [2000], Lunine [1999], and 
Yung and DeMore [1999]). 

Once we know the abundances of key gases, we 
can then make informed speculations on the likeli- 
hood that life exists on the planet. However even 
with the restriction to terrestrial- type planets, we still 
have a large range of possible types of atmospheres. In 
the following sections, we begin to explore the range 
of possibilities. We start with a description of our 
method for calculating visible and infrared spectra for 
the case of the present Earth. We extend this to in- 
clude Neoproterozoic icehouse and hothouse Earths, a 
methane-rich Earth, and the present Mars and Venus. 
We conclude with a brief discussion of how the pres- 
ence of life might be inferred from spectra such as 
these. 
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Figure 2. Solar system blackbody thermal emission 
spectra and reflected light spectra at 10 pc, for the 
Sun, Jupiter (J), Earth (E), Venus (V), Mars (M), 
and zodiacal dust (Z). For curve Z a telescopic field 
of view of 0.010 arcsec diameter centered at 0.1 arcsec 
from the star is assumed. For the Earth, present at- 
mospheric abundances are used to calculate a line by 
line spectrum for the entire 4 decades in wavelength, 
for a cloud-free atmosphere. (Note 1 Jansky is 10~ 26 
watt m -2 Hz -1 .) 

4. PLANETARY SPECTRA: 
PRESENT EARTH 

A graphical overview of the exoplanet detection is- 
sue is shown in Figure ||, where we plot the flux den- 
sity of a model solar system as it would be seen from a 
distance of 10 pc, the median distance for the nearest 
450 or so stars in our Galaxy. 

Blackbody spectra of the Sun, Jupiter, and the 3 
terrestrial planets are shown in Figure |[ for the ef- 
fective temperatures of these bodies, except as fol- 
lows: the spectrum of Jupiter includes its internal 
heat source contribution; the spectrum of Mars is an 
average of the day and night side spectra at differ- 
ent temperatures; and for Earth the average ground 
temperature is used, not the effective temperature. 
(N.B., the effective temperature of an object is the 
temperature of a blackbody which has the same area 
and total radiated power as does the object.) 



The zodiacal dust cloud is modeled as a face-on, 
smooth, optically thin blackbody emitter with optical 
depth varying as r~ ' 39 and temperature varying as 
r -° A2 where r is distance from the central star [Reach, 
1995]. 

The reflected light from Jupiter and the terrestrial 
planets is approximated by a scaled version of the so- 
lar spectrum, proportional to each planet's average 
visual albedo and area, reduced by a factor of 0.26 to 
approximate the brightness at quadrature, when the 
observer sees only one-half of the disk illuminated. 
The zodiacal dust reflection spectrum is modeled with 
the same density distribution as for the thermal emis- 
sion but with an albedo selected to agree with visual 
observations. 

The spectral line component of the Earth's spec- 
trum in Figure || is calculated separately for the ther- 
mal emission and reflection cases, for a clear atmo- 
sphere, and the results combined, as described next. 

4.1. Spectral Computation Method 

Model Earth spectra are calculated with our SAO 
code originally developed to analyze balloon-borne 
far-infrared thermal emission spectra of the strato- 
sphere [e.g., Traub and Stier, 1976; Johnson et al, 
1995], extended to include visible reflection spectra. 
The spectral line data base includes the large AFGL 
compilation [Rothman et al., 1998] plus improvements 
from pre-release AFGL material and our own sources. 
In a few cases laboratory cross section spectra are 
available but spectroscopic analysis is not, so here 
we use an empirical pseudo-line band shape. The 
far wings of pressure-broadened lines can be non- 
Lorentzian at around 1000 times the line width and 
beyond, so in some cases (H 2 0, CO2, N 2 ) we re- 
place the far wings of the line-by-line calculation with 
measured continua data in these regions. Dust and 
Rayleigh scattering are approximated by empirical 
wavelength power laws and contribute significantly 
only in the visible blue range. Model atmospheres 
from to 100 km altitude are constructed from stan- 
dard models discretised to appropriate layers, and 
additional radiative transfer methods used to ensure 
that line cores and optically thick layers are accu- 
rately represented. 

Radiative transfer from layer to layer is explicitly 
calculated using the average absorption and emission 
properties of each layer; scattering as a source is ne- 
glected. Integration from the spherical Earth atmo- 
sphere is approximated to a few percent accuracy by 
a single-point calculation at a zenith angle of 60 de- 
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grees, so the effective air mass is 2 in the infrared 
(outgoing emission) and 4 in the visible (2 for incom- 
ing sunlight, plus 2 for outgoing reflected light). 

Cloud effects are beyond the scope of this chapter, 
and they are not included in the calculations shown 
here, but they can be represented by inserting con- 
tinuum absorbing/emitting layers at appropriate al- 
titudes; broken clouds can also be represented by a 
weighted sum of spectra using different cloud layers. 
In general, the effect of clouds is to dilute the strength 
of line features in the visible, and to dilute, but in ex- 
treme cases cause absorption lines to appear as emis- 
sion lines, in the infrared. 

4.2. Thermal Emission Spectrum 

The dominant features of the Earth's thermal emis- 
sion spectrum are illustrated in Figure |L where the 
blackbody flux and composite spectrum (Jy/sr) are 
shown in the top left panel, and the other panels show 
the relative intensities of the major infrared molecu- 
lar species (H 2 0, O3, CH4, CO2, N 2 0) as well as 
minor contributors (H 2 S, S0 2 , NH 3 , SF 6 , CFC-11, 
CFC-12). The composite spectrum is calculated for 
the present abundances of each species, but the in- 
dividual species spectra are calculated for increased 
or decreased abundances, with the expected vertical 
mixing ratio profiles scaled so as to show the absorp- 
tion spectrum minima at an optical depth of about 
unity. 

The H 2 panel shows the far-infrared rotational 
band and the mid-infrared vibrational band, calcu- 
lated for an abundance of about 0.1 times the satu- 
rated value in the lower troposphere. Even at this re- 
duced concentration, the water lines are quite strong, 
but also rather diffuse, without any well-defined com- 
pact spectral features in the thermal infrared. 

Carbon dioxide, on the other hand, even at 5% of 
present abundance, shows a very strong 15 /zm (667 
cm -1 ) band in the infrared, the depth of which is 
limited not by abundance, but rather by the thermal 
structure of the Earth's atmosphere, such that the 
minimum brightness corresponds to the blackbody 
strength at the altitude at which the band core opti- 
cal depth reaches values on the order of unity. Note 
the small spike in the center which is generated at an 
altitude of about 30 km by the temperature inversion 
in the stratosphere, causing the very strong core to 
appear in emission against the lower-lying cooler lay- 
ers from which the near-wing emission emanates at 
an altitude of about 20 km. 
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Figure 3. A calculated thermal emission spectrum 
of the present Earth is shown in the top left panel. 
The other panels show spectra of individual species, 
with their mixing ratio profiles scaled up or down so 
as to generate a maximum feature depth of roughly 
one-half. The calculations are performed at very high 
resolution and subsequently smoothed to 1 cm -1 for 
display. The panels demonstrate the concentrations 
needed for each species to contribute significant opac- 
ities on an Earth-like planet. 

Ozone is the third most prominent infrared fea- 
ture with its strong 9 /xm (1100 cm -1 ) band, shown 
here at its natural abundance. The ozone feature 
is almost entirely due to the stratospheric O3 layer, 
though there is a trace of tropospheric O3. 

Methane and nitrous oxide are both shown at 
twice natural abundance, and have significant fea- 
tures nearly overlapping in the 7 /jm (1400 cm -1 ) 
region, also lying in the red wing of the 6 jum (1600 
cm -1 ) water band, and therefore not readily separa- 
ble, but nevertheless in principle measureable. The 
combined effect of the current abundance of CH4 and 
N 2 is seen in the composite panel, where the spec- 
trum shows a rather sharp decrease going from about 
1200 to 1300 cm -1 , compared to the gradual decrease 
in this region due to H 2 alone, as shown in the 
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water-only panel. 

The remaining species, at present Earth abun- 
dances, are not expected to be easily detectable on an 
exoplanet, due to the weakness of their bands. The 
panels show H2S at 10,000 times natural abundance, 
S0 2 at 1000 times, NH 3 at 1000 times, and the an- 
thropogenic gases SF 6 , CFC-11, and CFC-12 each at 
100 times current abundance. 

4.3. Reflection Spectrum 

Dominant features in the present Earth's reflec- 
tion spectrum are shown in Figure |^, where in this 
case the panels give the reflected intensity normal- 
ized to the incident solar intensity, smoothed from 
the original high-resolution calculation to a plotted 
resolution of 100. The upper left panel is a composite 
of the 5 contributing species H 2 0, O2, O3, CH4, and 
CO2, for present abundances. Note that the spectral 
range goes from the near infrared (2000 cm -1 = 5 
/mi) to the near ultraviolet (33000 cnr 1 = 0.30 /jm), 
with a continuum normalized everywhere to unity. 
We have ignored the thermal emission contribution 
at long wavelengths as well as the Rayleigh and dust 
scattering components which will show up mostly at 
short wavelengths. 

The H 2 panel, for the present abundance of wa- 
ter, shows a series of absorption bands spanning the 
middle part of the visible spectrum and increasing in 
strength toward the near infrared. The strengths of 
these vibrational bands are essentially independent 
of temperature, but will increase in proportion to 
the abundance of water and the square root of air 
pressure, however since the lines are relatively satu- 
rated, the average band depth will only increase as the 
square root of band strength. The net result is that 
these bands should be good indicators of the presence 
of water over a large dynamic range of conditions, al- 
though this same property makes them less useful as 
quantitative indicators of water mixing ratio, unless 
we also have independent knowledge of temperature 
and pressure. 

The strongest 2 band is the Fraunhofer A-band 
at 0.76 fim (13000 cm -1 ). This band too is satu- 
rated, and will still be relatively strong for signifi- 
cantly smaller mixing ratios than the present Earth's. 
It will therefore be an excellent indicator of the pres- 
ence of oxygen (see Chapter VI. 2). 

The O3 molecule has two broad features of note, 
the extremely strong Huggins band which produces 
the ultraviolet absorption shown here at about 0.33 
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Figure 4. The reflectivity of the Earth with present 
atmospheric abundances is shown in the top left 
panel, normalized to unity. Clouds, aerosols, and 
Rayleigh scattering are ignored in this example. The 
other panels show reflection spectra for the cases 
where only a single species is present in terrestrial 
abundances. For CH4 and CO2 we also show spectra 
for greatly enhanced abundances as discussed in the 
text. 

/im (30000 cm" 1 ) and shorter, and the Chappuis band 
which shows up as a broad triangular dip in the mid- 
dle of the visible spectrum from about 0.45 to 0.74 
/tm (22000-13000 cm" 1 ). Ozone in the stratosphere 
is produced from O2 molecules, and its abundance 
is a non-linear function of the O2 abundance, such 
that even a small amount of O2 can produce a rel- 
atively large amount of O3 [Kasting and Donahue, 
1980]; however in absolute terms the modeled column 
abundance of O2 is nevertheless large compared to 
that of 3 by a factor of 20,000 to 500,000. 

Methane at present terrestrial abundance (1.65 
ppmv) has no significant visible absorption features, 
but at high abundance (0.1%) it has strong visible 
bands at 0.9 and 1.0 /mi (11000 and 10000 cm" 1 ). 

Carbon dioxide has negligible visible features at 
present abundance, but in a high-C02 atmosphere 
(10%) it has a significant band at 1.2 /im (8000 cm" 1 ) 
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and even stronger ones at longer wavelengths. 

5. PLANETARY SPECTRA: 
NEOPROTEROZOIC EARTH 

The present status of our secure knowledge of 
paleo-Earth atmospheres is easy to review, because 
little of our knowledge is secure, although in recent 
years the situation has been improving. In overview, 
the evolution of CO2 has long been believed to have 
decreased from a high level of roughly 1 bar at about 
4.5 Ga (where Ga represents 1 x 10 9 years ago) to 
0.00035 bar at present, however recent evidence sug- 
gests that major oscillations occured around 0.5-0.8 
Ga, and perhaps at other glaciations, The abundance 
of O2 is believed to have been less than 0.001 bar until 
about 2 Ga when it rapidly began to increase toward 
its 0.2 bar present value, and the time of transition 
is roughly coincident with the beginning of abundant 
phytoplankton on Earth. 

The equilibrium temperature of the Earth with its 
present albedo, present solar flux, and no greenhouse 
gases is about 246 K, below the freezing point of wa- 
ter; this corresponds to the physical temperature at 
a level in the stratosphere at which the Earth effec- 
tively radiates. The Earth is rescued from freezing 
by the greenhouse effect, which can be envisioned as 
being driven at temperatures below the freezing point 
of water by CO2 and aided at higher temperatures by 
evaporated water vapor. The combination of present 
levels of CO2 and H2O is sufficient to warm the sur- 
face to 290 K, about the current average surface tem- 
perature, which is of course sufficient for liquid water. 

The early Earth was illuminated by a weaker Sun, 
about 0.71 times the present luminosity, which with- 
out CO2 would have resulted in an even cooler Earth, 
about 246(0.71) 1/4 = 226 K. The argument for a large 
amount of early CO2 is simply that we believe that 
liquid water, with only intermittent glaciation, was 
present during most of Earth's past [cf. e.g., Lunine 
1999, Sec. 11.10 and 19.4], and this requires a large 
greenhouse effect, equivalent to about 0.2 bar or more 
CO2. Also, since there is the equivalent of about 50 
bar of CO2 deposited in crustal rocks, the pressure 
may have been higher. 

The level of CO2 probably did not fall steadily 
since 4.5 Ga. We know that periods of major glacia- 
tion did occur, and that almost certainly these could 
not have begun unless the CO2 abundance had first 
dropped significantly. Recently there has been a ma- 
jor advance in this area with the clear-cut identifi- 



cation of 3 major glaciation cycles alternating with 
warm tropical conditions at about 0.5-0.8 Ga [Hoff- 
man et al., 1998]. This period is roughly the end 
of the Proterozoic era (from 4.5 to 0.7 Ga), and the 
start of the Phanerozoic era (from 0.7 to 0.0 Ga), 
so the sudden glaciations are said to have produced 
a Neoproterozoic "snowball" Earth. The CO2 abun- 
dances probably oscillated between roughly 0.2 bar 
and 0.0001 bar. The evidence for oscillation is un- 
ambiguously seen in rock layers in Namibia: repeated 
pairings of rounded boulders (from glaciers) topped 
by thick layers of carbonate rock (precipitated from 
seawater). 

A major mystery is why these oscillations have not 
occurred over most of Earth's history, or more point- 
edly, why the oscillations started when they did and 
stopped after a few cycles. It is possible that the other 
major glaciations were accompanied by such oscilla- 
tions, but this has not yet been established, or even 
discussed, to our knowledge. 

One might expect that the icehouse state would be 
triggered by a sudden drop in CO2, caused perhaps by 
a rapid uptake of CO2 in the ocean, for some reason, 
and that the hothouse state would be triggered by the 
accumulation of volcanic CO2, a direct result of plate 
tectonics, which at present rates would take only a 
geologically short interval of about 10 million years 
to reach greenhouse levels. However the details of the 
oscillations are as yet not well understood. 

The implication of the Neoproterozoic oscillations 
for extrasolar planet searching is clear: we might en- 
counter such a condition on an extrasolar planet, and 
we must be prepared to anticipate its spectral signa- 
ture. We show in Figure |^ our preliminary calculation 
of the infrared spectra of a Neoproterozoic hothouse 
and icehouse Earth. To simulate the icehouse state, 
we use the present Arctic temperature profile, and O2 
at 1% PAL (present atmospheric level); the O3, H2O, 
N2O, and CH4 are modified [Kasting, 1980] for a 1% 
O2 level, and the CO2 is set to 100 ppm; about 1 bar 
of N2 is also assumed to be present. To simulate the 
hothouse state, we used a present tropical tempera- 
ture profile, O2 at 10% PAL, and corresponding O3, 
H 2 0, N 2 0, and CH 4 ; also C0 2 is set to 120,000 ppm. 

A substantial literature is building on the snow- 
ball Earth concept, for example climate simulations 
[Jenkins and Smith, 1999; Hyde et al., 2000] and re- 
alizations that life would not be totally eliminated by 
ice cover [Runnegar, 2000]. 
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found at another epoch, 0.183 Ga [Hesselbo et al, 
2000] . This event is associated with high surface tem- 
peratures and significant mass extinction. Yet other 
such events at 0.090 and 0.120 Ga have been similarly 
reported [Kerr, 2000]. Interestingly, each of these 
events, now 5 in total, occurs at about the same time 
that a large volume of volcanic outflow occured, and 
a corresponding mass extinction occured, suggesting 
that perhaps the volcanoes triggered the methane and 
this in turn triggered further warming which resulted 
in large-scale loss of species. Large amounts of CH4, 
such as in these bursts, will produce measurable ab- 
sorption features in both the infrared emission and 
near-infrared absorption spectra of the Earth. 

7. PLANETARY SPECTRA: 
PRESENT VENUS AND MARS 



Figure 5. A calculated thermal infrared emission 
spectrum for a Neoproterozoic hothouse condition 
(left), and icehouse (right). Note the huge difference 
between these; for example, the 15 /urn CO2 feature 
is an apparent emission feature in the hothouse spec- 
trum, but it is an apparent absorption feature in the 
icehouse spectrum. The apparent emission feature is 
a result of a warm inversion layer in the stratosphere 
(due to O3 heating) combined with a CO2 mixing ra- 
tio which gives an optical depth of about unity at the 
top of this layer for wavelengths in the core of the 
C0 2 band. 

6. PLANETARY SPECTRA: 
METHANE 

Besides CO2, methane (CH4) could strongly influ- 
ence the surface temperature of a planet. Methane 
bursts may have punctuated the Earth's history. For 
example, carbon and oxygen isotope fractionation in 
marine carbonate deposits strongly imply that about 
350 times the present level of CH4 was injected into 
the Earth's atmosphere, during a span of less than 
1000 years, about 0.055 Ga [Bains et al, 1999]. The 
source of the CH4 may well have been methane hy- 
drate, which is found on the sea floor in abundance 
even today. That this time of injection coincided 
closely with a period of warming (the Late Paleocene 
thermal maximum) strongly suggests a cause and ef- 
fect relationship between these phenomena, but it is 
not clear which is cause and which is effect. 

Further evidence for a methane burst has been 



We have calculated infrared spectra of Venus and 
Mars, using current models of atmospheric abun- 
dances and temperature profiles, with the added fac- 
tor that the opacities of air-borne dust and aerosol 
droplets can be included in order to reproduce ob- 
served silicate and ice features. Clearly we have much 
less information on the geologic history of these plan- 
ets than we do for Earth, but spectral evidence of 
isotopic fractionation on both planets has led to in- 
formed speculation about the past abundance of wa- 
ter, for example. 

We show in the right panel of Figure || our cal- 
culated Mars spectrum using a recent model atmo- 
sphere [Y. Yung, 1999, private communication]. The 
CO2 and H 2 features in this spectrum are a good 
match to many observed Mars spectra. We have re- 
peated these calculations with the inclusion of thin 
ice clouds and silicate dust clouds, and these pro- 
vide a reasonably good simulation of certain other 
observed Mars spectra. The calculations for Venus 
are shown in the left panel of Figure [|, based on the 
model atmosphere of Bullock and Grinspoon [2001]. 
This spectrum is dominated by CO2, and does not 
include any aerosols, so the fact that the very bright 
thermal emission from the ground is not seen at the 
top of the atmosphere is entirely due to the opacity 
of the wings of the strong CO2 bands. 

8. CAN WE DETECT LIFE? 

The big question, can we detect life on an ex- 
trasolar planet, was addressed in a general sense by 
Lovelock, [1965] who advocated searching for signs of 
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Figure 6. Calculated spectra for Venus and Mars. 
Both spectra are dominated by CO2. No ice or dust 
features are included here, although these features are 
seen occasionally on Mars, and we have modeled them 
successfully. 

chemical non-equilibrium, such as the simultaneous 
presence of reducing and oxidizing gases. In princi- 
ple, if one knew the abundances of all species present 
in an atmosphere, and the boundary conditions (in- 
cident spectrum from star, surface composition, ro- 
tation rate, etc.), then one could calculate the state 
of photochemical and dynamical equilibrium, and its 
likely fluctuations, and compare this with the ob- 
served state, to see if there are significant differences 
which might be interpreted as signs of life, a term 
not included in the equilibrium calculation [cf. Nis- 
bet and Sleep, 2001; Brack, 1998]. However our expe- 
rience with the atmosphere of the present Earth tells 
us that the observed state is frequently not predicted 
by theory, so this approach may not be foolproof 

From a spectroscopic point of view, it is natu- 
ral to think of searching for life in terms of measur- 
ing the abundances of atmospheric constituents, and 
comparing these abundances with the results of cal- 
culations of thcrmochcmical and photochemical equi- 
librium. In particular, the simultaneous presence of 
significant amounts of oxidized and reduced species, 
such as H2O, CO2, N2O, and CH4, or the presence 
of large amounts of O2 or O3, would both be indi- 
cators of life. Indeed both conditions are found on 
Earth, where there is no known means of producing 



anything approaching a 21% O2 atmosphere except 
by photosynthesis, and the main sources of N2O and 
CH4 are biological. Earth's O3 is produced in rela- 
tively large amounts in the stratosphere by photolysis 
and recombination of O2, making it a good indicator 
of the presence of O2 . As signs of life in themselves, 
H 2 and CO2 are secondary in importance, because 
although they are raw materials for life, they are not 
unambiguous indicators of its presence. Farther down 
the chain, CH4 is a life product on Earth, but else- 
where it is also a ubiquitous primordial species, and 
it is difficult to detect spectroscopically unless it is 
quite abundant. Likewise, N 2 is interesting because 
it is produced in abundance by life, and only in trace 
amounts by natural processes, but it can only be de- 
tected in a region which is strongly overlapped by 
CH4 and H2O, and so is an unlikely prime target. 

As counterexamples, we know that small amounts 
of O2 and O3 are readily produced by non-biological 
means. For example, a trace amount of O2 is seen 
on Venus, consistent with photochemical production 
from CO2; on Mars both O2 and O3 are measured, 
but the amounts are in agreement with calculations 
of photochemical production from CO2 and loss due 
to CO2 and H 2 and products [e.g., Yung and De- 
More, 1999]. In two other examples, Europa has 
about 10~ n bar of O2 in its atmosphere (e.g., Chap- 
ter III. 3), and Ganymede may have O2 trapped in 
surface ice, but in both cases the observed signatures 
are consistent with energetic particle bombardment of 
a water-ice surface, and therefore they do not require 
a biological source. 

An intriguing theory of the origin of life on Earth 
[Wachtershauser, 2000] has found experimental back- 
ing in new laboratory experiments [Cody et ai, 2000] 
which show that iron sulfide at elevated pressure and 
temperature can facilitate the natural generation of 
organomctallic compounds by an autocatalytic pro- 
cess, in which the products of reaction catalyze the 
next cycle of reaction, making the process likely to be 
self-sustaining. This work is of interest because iron 
and sulfur atoms are at the center of many enzymes, 
and the process provides a simple explanation of their 
presence in living cells. The work is also of interest 
because it suggests that the origin of life on Earth 
may have been deep underground, where conditions 
are right for the reactions to take place, and out of 
the possibly less-friendly conditions on the surface of 
the early Earth. 

The Neoproterozoic oscillations discussed above 
may also be critically important to the development 
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of life on Earth, because, as the eponymic oscillations 
suggest, this era coincided with the explosive radia- 
tion of life forms on Earth, from simple algae types to 
the dramatically different types of life forms we know 
today. It is thought that the huge stresses, due to 
oscillations of temperature, nutrients, and available 
sunlight, imposed on simple life forms forced them to 
develop new adaptations in order to survive in their 
new environments, in what may have been the great- 
est Darwinian experiment ever. Thus, if we find ex- 
trasolar planets in the throes of such oscillations, we 
may expect that a similar life-radiation may be oc- 
curing on that planet, although we will have to wait 
another several tens or hundreds of million years be- 
fore we can see how the experiment turns out. 
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